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Abstract: 

We propose an asymptotic expansion formula for matrix integrals, including oscil- 
latory terms (derivatives of theta-functions) to all orders. This formula is heuristically 
derived from the analogy between matrix integrals, and formal matrix models (combi- 
natorics of discrete surfaces), after summing over filling fractions. The whole oscillatory 
series can also be resummed into a single theta function. We also remark that the co- 
efficients of the theta derivatives, are the same as those which appear in holomorphic 
anomaly equations in string theory, i.e. they are related to degeneracies of Riemann 
surfaces. Moreover, the expansion presented here, happens to be independent of the 
choice of a background filling fraction. 
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1 Introduction 



Convergent matrix integrals of the form 




dM e 



(1.1) 



are very usefull in many areas of physics (statistical physics, mesoscopic physics, quan- 
tum chaos,...) and in mathematics (probabilities, orthogonal polynomials,...) [TT]. 
People are mostly interested in their asymptotic behavior in the large n limit (and 
n/N ~ constant). 

There is another form of matrix integrals, called formal-matrix integrals, which 
come from combinatorics (2d quantum gravity for physicists [TDl [151 [IB])- They are 
generating functions for counting discrete surfaces (also called "maps") of given topol- 
ogy. Formal matrix integrals are only asymptotic series, they are not convergent in 
general, and almost by definition, they always have a large n expansion (see All 
the terms in their large n expansion are known pTl [22] , and are deeply related to alge- 
braic geometry and integrable systems. They have many applications to combinatorics, 
and string theory in physics [3ll [36] . 

In this article, we use the analogy between the two types of matrix integrals, and 
generalizing the method of [D] , we propose an asymptotic formula for convergent matrix 
integrals, including oscillations to all orders: 



where B is a theta function, i.e. a periodic function, this is why we call G and its 
derivatives "oscillatory terms". 

Then we observe that the series containing the oscillatory terms can be resummed 
into a single theta function: 



We also observe that the coefficients in front of the derivatives of in eq ll.2l are 
the same which appear in the so-called "holomorphic anomaly equations" discovered 
in the context of topological string theory |8j. In other words they are related to the 
combinatorics of degeneracies of Riemann surfaces. 



Z 




(1.2) 




(1.3) 



2 



Finally, we observe, that although we define each term of the expansion after choos- 
ing a reference filling fraction e*, the partition function is in fact independent of that 
choice. This is related to the so-called background independence problem in string 
theory, first observed by Witten [55] . 

For the 1-hermitian matrix model (with real potential), the first term of this asymp- 
totic expansion 

^^e^«^'"'^^« e (1.4) 

was derived rigorously by Deift& co [1A\ using Riemann-Hilbert methods, and their 
method proved the existence of a whole oscillatory series containing derivatives of the 
Theta-f unction. The same result was also obtained by heuristic physicists methods by 
[9]. Here, we generalize the method of [9] and we give the exact coefficient of the whole 
series. 

Also, in case where the genus of the Theta function is zero, there is no oscillatory 
term, and one finds the so-called topological expansion Z ~ e^f^^ ^^^3^ which is 
well known to coincide (in the sense of asymptotic formal series) with the generating 
function for enumerating discrete surfaces [10] . In this genus zero case, the asymptotics 
of the convergent matrix integral were derived by several methods and several authors 
\W\ [50] . The coefficient of the expansion are of course the symplectic invariants of [22] • 

For other convergent matrix models, for instance the 2-matrix model, such expan- 
sions were conjectured many times [231 [21|, but never proved. Here, we don't prove it 
either. We merely give all the coefficients in the formula to prove, and we explain their 
heuristic origin. 

As we said above, the heuristic origin of the formulae presented in this article, is 
just the analogy between formal and convergent matrix models. 

Outline: 

• In the first section, we define the convergent matrix model on generalized paths, 
and write it as a sum over filling fractions. 

• In the second section, we consider the formal matrix model. 

• In the 3rd section we perform the sum over filling fractions, and we get 0- 
functions. 

• In the 4rth section we discuss the link with degenerate Riemann surfaces and 
holomorphic anomaly equations. 

• In the 5th section we discuss the background independence problem. 

• Section 6 is the conclusion. 
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1.1 Introductory example: 1 matrix model 
1.1.1 Paths and homology basis 

Consider a polynomial potential V{x), of degree d + 1 = deg\^, with complex coeffi- 
cients. There are many different paths 7 such that the integral 

dxe-^^") (1.5) 

7 

is absolutely convergent. These are the paths which go to cxd in a sector where Re^ > 0, 
or more precisely, the paths which connect two such sectors (see [5] for a discussion on 
that, or |3]. Those considerations can be easily extended to any V such that V is a 
rational fraction). 

Example: quartic potential V{x) = x'^, we have deg\^ = 4, i.e. there are d = 3 





independent paths, for example we choose: / I ^ . In this example, 
we have M = 72 + 73. 

In fact, there are d = degV homologically independent such paths. Let us choose 
a basis: 

7i,---,7rf (1-6) 

This means any (unbounded) path on which the integral dx e~^*^^^ is well defined, 
can be decomposed on the basis: 

d 

7 = ^Ci'yi (1.7) 

By definition: 



1=1 



! dxe-^(^) = ^Q I rfxe-^(^) (1.^ 
•Ji i=i J■y^ 



In this definition, the Cj's can be arbitrary complex numbers, they don't need to be 
integers. 

However, if 7 is itself a path, the Cj's can take only the values +1,-1, or 0. 
If the Cj's are not integers, we say that 7 = Yli Cj7i is a generalized path. 
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1.1.2 Matrix model on a generalized path 

Let 7 be a generalized path. We define the set of Normal matrices on 7: 



//„(7) = {M = C/diag(xi, ...,Xr,)U^ / U e U{n) , Vi e 7} (1.9) 
equipped with the measure: 

dM = A(x)2 dU dxi ... dxn , A{x) = Yl{xj - Xi) (1.10) 

where dU is the Haar measure on U{n), and A(x) is the Vandermonde determinant, 
and dx is the curviline measure along the path (if 7 = a;(s) , s e M, is a parametrization 
of the path we have dx — x'{s) ds). 

Remcirk: Hn(^) — is the set of hermitian matrices, with the usual U{ri) 
invariant measure. 

Remcirk: In general Hn{'y) is not a group, for instance the sum of two matrices in 
Hni'f) is not in i?n(7), and the product by a scalar is not either. Also, the "measure" 
dM is not positive, in fact it is complex. 

The matrix integral on i?„(7) is defined as follows: 

Z(7) = i / dMe-^^^(^) = l / dx,...dxnA{x)^ Ue-''''^-^^ (1.11) 
n\ Jh„{j) ril J^n Y 

or in other words: 

ni+...+nd=n 

where we have defined: 

Z{n,/N, na/N) = — ^ / dx,... dx^ A{x)' TT e-^^^^^) (1.13) 

nil...ndl J^^^x...x^2'' t 

1.1.3 Assumption: topological expansion 

First, let us assume that only ^ + 1 < d of the q's are non-vanishing. We write: 



yi^l,...,g, a^e^''"'' , Cg+i^l , yi^g + 2,...,d, q = 



(1.14) 

If 7 is a path, the q's take the values ±1, and thus Ui — or 1/2. The vector {ui, . . . ,1^) 
is going to be considered a charcateristic in a genus g Jacobian. Also, up to reverting 
the orientation of 7^, we can always assume that if 7 is a path, 

7 = path ^ Vi = 1, . . . 1, Q = 1 , ^ i/ = (1-15) 
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Hypothesis: 

Our working hypothesis is that the basis paths 71, . . . ,73+1 have been chosen so 
that each Z{ni/N, . . . ,nd/N) admits a large topological expansion: 

00 

In (Z(ei, . . . , e,_i)) ~ F{e) = N'-'^F^'^\e) (1.16) 

h=0 

It is conjectured that given a (generic) potential V, and a generalized path 7, such 
a "good" basis always exists (may be not unique). In fact, for the 1-matrix model 
with polynomial potential, this can be proved a posteriori from the asymptotics of M. 
Bertola [6l[7]. But for more general cases, it is only a conjecture, for instance for the 
2-matrix model. 

Now, let us explain where this hypothesis comes from, and what heuristic arguments 
support it. 

1.1.4 Loop equations and Virasoro constraints 

It is well known that any integral defined in eq J1.13l satisfies an infinite set of lin- 
ear equations, sometimes called "loop equations" [15], or Virasoro constraints, or 
Schwinger-Dyson equations, or Euler-Lagrange equations, and which just come from 
integration by parts: 

VA:>-1, Vk.Z = (1.17) 

They satisfy Virasoro algebra: 

[Vk,V,] = (k - j)Vk+j (1.19) 

Remark: It is important to notice that, since integration by parts is independent 
of the integration paths (as long as there is no boundary term), both Z{'j) and any 
Z{ni/N, . . . ,nd/N), Vra,, satisfy the same set of loop equations. 

1.1.5 Formal matrix models and combinatorics of maps 

Formal matrix integrals are defined as formal generating functions for enumerating 
discrete surfaces (also called "maps", i.e. topological graphs embedded on a Riemann 
surface, such that each face is a disc) of given topology. Basically, Fg is the generating 
function for counting "maps" of genus g. The generating series: 

00 

lnZformai = $^iV'"'^F3 (1-20) 

5=0 
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needs not be convergent, and in fact it is never convergent if the weights for "maps" 
are positive. It is merely a formal series, whose only role is to encode the -Fg's. 

The formal matrix integrals satisfy the same loop equations, i.e. Virasoro con- 
straints as Z{'~f) and Z{ni/N, . . . ,na/N) (see [U]). In the context of combinatorics of 
maps, loop equations are known as Tutte's equations [3Z], and were first obtained by 
counting "maps" recursively (removing one edge at each step). 

The Fg's of formal matrix integrals have all been computed: Fq has been known 
for a long time, then Fi [12], and all the F^'s with g > 2 were computed recently in 
[201 [22]. 

In fact, it was proved in [201 [22], that any solution of loop equations which has a 
topological large N expansion of the form: 

oo 

InZ = J2N^~^^Fg (1.21) 

9=0 

can be obtained by the symplectic invariants of ^2], i-e. they are encoded by a spectral 
curve. 

1.1.6 Spectral curves 

Both the convergent matrix integral, and the formal matrix integral are associated to 
an (algebraic) spectral curve of the form: 

y^ = V"W^-l/lV^^M^IM\ (1.22) 

^ ^ ' n\ x-M I ^ ' 

• For the convergent matrix integral Z defined in eq jl.lll the average < . > is 
taken with respect to the measure dM e-^Tri/(M)_ rpj^g notion of a spectral curve, 
comes from the orthogonal polynomials method of Dyson-Mehta [35J, combined with 
the theory of integrable systems [2]. The orthogonal polynomials satisfy an integrable 
differential equation of the form = T>[x) ip, where V{x) is a 2 x 2 matrix with 
polynomial coefficients, and the spectral curve is by definition the set of eigenvalues of 
V (Jimbo-Miwa-Ueno [32]), i.e.: 

y2 = lTrI?(a;)2 (1.23) 

It was proved |4j that: 

2 ^ ' ^ ' N\ x-M / ^ ' 
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• For the formal matrix model, and more generally, for an arbitrary solution of the 
Virasoro constraints which has a topological expansion, the average < . > has a formal 
meaning, and can be defined from the Virasoro generators Vk- It is not the purpose 
of this article to explain where it comes from (see [IHl US]), and the spectral curve is 
the algebraic equation satisfied by the "disc amplitude", i.e. generating function for 
counting planar "maps" with one boundary (i.e. having the topology of a discs), and 
it can be proved that it satisfies: 

= V'{x)^ -AP{x) (1.25) 

where P{x) is a polynomial of degree d — 1 = deg\^", and with the same leading 
coefficient as V. 

d-2 

P{x) = {d + l)U+ix''-' + Y,Pkx'' (1.26) 

fc=0 

The coefficients Pk, are the conserved quantities in the context of integrable systems 
[2], whereas the coefficients of V are called the Casimirs. The coefficients Pk are in 
1-1 correspondance with the so-called "action variables": 

^i = W~ f ' i = 1,. . . ,d-l (1.27) 

Here in the random matrix context, the ej's are called filling fractions. 

1.1.7 Symplectic invariants 

In [2^, it was proved, that given a spectral curve 

S{x,y) = (1.28) 

(here S{x,y) = — V'^xY + 4P(x), i.e. in other words, given a potential V{x) = 
X^fcti tkx'^ ^iid a polynomial P{x) = {d+ 1) t^+i x'^~^ + J2k^ oi in other words 

given V and the filling fractions e^'s), it is possible to define an infinite sequence: 

Fg{£) , ^ = 0,...,oo (1.29) 

such that: 

oo 

t{£) = exp J2 N^~^'Fg{£) (1.30) 

9=0 

is a solution of loop equations. 

The Fg{£) were constructed in |22j for any spectral curve £{x,y) = 0, and they 
have many interesting properties, for instance they are invariant under symplectic 
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deformations of the spectral curve, and t{S) is the r-function of an integrable hierar- 
chy. Their modular properties were also studied in |22j and further clarified in [19] , and 
they happen to be deeply related to the so-called Holomorphic anomaly equation 
first found in string theory [SI |T], and which relate the non- holomorphic part of the 
generating function for counting Riemann-surfaces to the contribution of degenerate 
Riemann surfaces (nodal surfaces). This will play a role below. 

Also, in [22], were defined the correlators: 

W^'\z,,...,z,) , g = 0,...,oo , k = 0,...,oo , { W^'^ = F, ) (1.31) 

which are multilinear symmetric meromorphic differential forms on the spectral curve. 
They also have many interesting properties, in particular they can be used to compute 
derivatives of the F^'s with any parameter on which S may depend. For instance 
derivatives with respect to filling fractions are: 

^ Wi'\z,, ...,Zk)= i ...,Zk, Zk+i) (1.32) 

(where r is the Riemann matrix of periods of the spectral curve, and Ai fl Bj = 6ij is 
a symplectic basis of non contractible cycles, see |27J [28] for algebraic geometry). 

1.1.8 Heuristic support to the conjecture 

The conjecture is supported by the following facts: 

• Both the convergent matrix integral Z{j^, . . . , 0, . . . , 0) defined in eq J1.13l 
and the formal matrix integral ^formal (ei, ■ ■ ■ ,^g) satisfy the same loop equations. 

• Since loop equations are linear, the space of solutions is a vector space. 

• For given V, both the convergent integral Z(^, . . . , 0, . . . , 0), and the formal 
^formai(ei5 • • • 5 ) are Specified by the same number of parameters, i.e. g parameters 
(indeed ni + . . . + n^+i = n, so that only g of them are independent). 

Those observations support the idea that there exists a good basis of the vector 
space of solutions, such that each basis function is at the same time formal and conver- 
gent, i.e. there exists a set of basis paths 7^, such that Z{j^, . . . , ^^^7^, 0, . . . , 0) admits 
a topological expansion. 

We do not prove this conjecture in this article, but we take it as an asumption. 
1.2 Generalization 2-Matrix model 

All this can be extended to a larger context, for instance the 2-matrix model, or the 
chain of matrices, or the matrix model coupled to an external field. 
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2 matrix model 

Consider 2 polynomial potentials Vi and V2, such that degVi = di + l,degV2 = 
d2 + 1. There are di x ^2 independent paths on C x C on which the following integral 
is absolutely convergent: 

did2 

dx dy e-^'^^^-^'^y^-^^y , 7 = J]Q7i (1.33) 
■y i=i 

where each 7^ is a product of a path in the x— plane and a path in the ?/— plane. 

Then, similarly to the 1-matrix case, we can also define a matrix integral on a 
generalized path (see [23j): 

Z(7) = / dM,dM2 ^-NTrmM,)+V2{M2)-KhM,) (^ 34) 

which satisfies: 

Z{l) = J2 ^1' • • • Z{n,/N, njN) (1.35) 

where we have defined: 

Z{ni/N, rid/N) = — / dxi Adyi ... dxn A dyn 

ni\...nd\ J^^w,„y,^^^<i 

A{x)A{y) JJe-^(^i("^)+^^(J^»)-"^«'') (1.36) 

i 

The 2-matrix model generalized integral satisfies loop equations (which form a 
W-algebra instead of Virasoro), which also come from integration by parts, and are 
independent of the path. In particular, each Z{ni/N, . . . ,nd/N) satisfies the same loop 
equations. 

There is also a formal 2-matrix model, which was introduced as a generating func- 
tion for bi-colored discrete surfaces, it was called the "Ising model on a random lattice" 
pi] . Almost by definition, the formal 2-matrix model has a topological expansion: 

In Z = J] Ar2-25^^ ^^37^ 

9 

The formal 2-matrix model satisfies the same loop equations as the convergent one, 
and the solution of loop equations was found in [211 [ISl [22], and it was found that the 
Fg^s are again the symplectic invariants of 

matrix model with external field 
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The same features also hold for the matrix models in an external field. The famous 
example is the Kontsevich integral [33j, also called "matrix Airy function": 



^Kontsovich 



dMe-^^^"^-*^^ (1.381 



whose topological expansion is the combinatorics generating function computing inter- 
section numbers. 

Summary 

In all cases, there is a convergent matrix model defined on generalized paths, and 
there is a formal matrix model which computes the combinatorics of some graphs. Both 
the convergent and formal model obey the same set of loop equations. 

The formal model has a topological expansion 

InZ = ^A^2-29^^ ^139^ 

a 

where the F^'s are the symplectic invariants of [22j|, computed for some algebraic spec- 
tral curve £{x, y) = 0. And in all cases the dimension of the homology basis of paths on 
which the integral is absolutely convergent, is the same as the genus 'g of the spectral 
curve, i.e. the number of filling fractions: 

7 = ^Ci7i 4^ ei = — (b ydx , i = 1, . . . ,g (1.40) 
i=i -^A 

In all those cases, the method we describe below should work. 



2 Formal matrix model 

Now, assume that Z{ei, . . . , ed-i) has a topological asymptotic expansion: 

oo 

In (Z(ei, . . . , e,_i)) = F(e) = ^'"''^/^(e) (2-41) 

h=0 

Each Fh must then be a solution of formal loop equations, and therefore it is given by 
the formulae of [22], and therefore each is analytical in the e^'s. 

Then, we choose arbitrarily a "prefered" filling fraction e*, and perform the Taylor 
expansion: 

F^(') = Y.^/t'^('-'*^' ' ^i'^ = ^(0 (2.42) 
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(k) 

Remark: We don't write the indices for readability, but is a tensor. For read- 
ability we write the formulae as if there were only one variable e, i.e. ^ = 1, but in fact 
we mean: 



k 

k=o ii,...,ik i=i --n ■ ■ ■ ^^ik 



oo ^ k 



(2.43) 

but for simplicity we shall write eq l2.42[ and we leave to the reader to restore the 
indices if needed. 

The derivatives of Fg, are given by eq jl.321 (see [22]): 

F,= <f Wt'\z„...,z,) (2.44) 



In particular, it is well known (see [22]), that 

F^= J) ydx (2.45) 
Jb 

and ^Fq = r is the Riemann matrix of periods (see [271 [25] for introduction to 
algebraic geometry) of the specral curve £, i.e. 

2^ IJI = = = £ ^'-''^ 

where (iMj is the normalized basis of holomorphic differentials [271 [28] on £^: 



duj = Sij (2.47) 



And thus we have (formally): 

Z{e) = Z(e*)e*''^'('"'*)^('"'*)e^*''^'^('"'*)^^^ 

— F^''^ . . . fS'^ it - e*)^'' (2.48) 

where we the sum carries only on > 1 and 2 — 2/ij — /j < for all i. 

One should notice that the exponential is now at most quadratic in the e's. 

3 Oscillations 

Now we are going to perform the sum of eq jl.121 

ni4-...+n5+i=n 
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where 

7 = 5^Q7i , Ci = e^'^''^ (3.50) 

i 

Since the filhng fractions = ^ take integer values (up to a factor), we have 
to perform a sum of exponentials of square of integers. Such sums are called theta 
functions. They play a key role in algebraic geometry. Let us recall a few properties 
[271 [28]. 

3.1 Theta functions 

We define the B-function: 

It clearly satisfies: 



96 _ 
dt du"^ 

It is related to the usual Jacobi-theta function: 



(3.52) 



ti f 

e{u, t) = ^^-^..,.(-^, -) e^^-^^^* (3.53) 

ZtTT m 

where (— A^e*, v) is called the characteristics. The Jacobi theta function with charac- 
teristics (a, h) is defined by: 

e^ i,{u, r) = ^ g2i7r{n+a)(«+fe) ^..(n+a)r{n+a) ^ q^^^^^ + 6 + ra, t) e*™™ ^^i-aiu+b) (3^54) 
n 

It takes a phase after translation along an integer lattice period n + rm: 

3.2 Convergent matrix model 

We thus have: 

Z(7) ~ J^e^-^'^ZformaKn/iV) 

n 

~ Y,cT ...of Z{n,/N,...,ng/N,Q,...,Q) (3.56) 

n 

The sum carries on integers n.i>Q and ^ • rzj = n. Therefore %+i = n — ^2^=1 '^j is not 
independent from the others. Another remark, is that in that sum we expect that only 
the vicinity of some extremal rij will dominate the sum, and that values of the n^s far 
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from the extremum should give an exponentially small contribution. That asumption 
allows to extend the sum to rii G Z. 

Then, we use the Taylor expansion of eq j2.48[ and we find (again we use tensorial 
notations): 



Z(7) ~ Z{e*) J2 



^2iTTY,i>^ini ^iTT{n-Nt*)T{n-Nt*) ^2inN({n~Nt*) 



k li>0 - , ^ " 



(3.57) 

where we recognize the B-function and its derivatives 



^(7)-^(.-)EE E 0^ 



(3.58) 

This formula is the main result presented in this article. 
For instance the first few terms in powers of are: 



/ 1 F'" 

z{i) - z{e*)(e + -{F;e' + -^e"'^ 



2 2 24 6 12 ' 



+ . 

(3.59) 

3.3 Resummation 

The expansion of formula J3.58I can be resummed into a single 0-function. We want to 
write it as: 

Z(j) = Z{e*)Q(u,t) (3.60) 

where 

oo oo 

iVF^ + ^A^i-^V^) , t = z7rr + ^Ar-2^tW (3.61) 



OO 

-TV T T-t/ 

U 

h=l h=l 

For instance, one easily finds the first orders: 



^(1) = + ^ V . ; 3 g2 

^ 6 e'(u(o),i7rr) ^ ^ 

,(1) _ , e----(^<°\^vrr) /G^ e'"\ (F^-Q^ /e(^) e---^\ 

2 24 e"(M(o),?7rr) 6 V ©" ©' / 72 V ©" ©'V 
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The Taylor expansion of eq J3.60l reads (and we use eq J3.52l) : 
Z(7) = Z(Oe(A'-Fo + ^«'" + ...,i»-r+ij(m + ...) 



^(^•) E E E 



m! n! 

m,n ki,...,km ji,---,jn 

fjm+2n 

^(fci) . . .^(fc^) tin) . . .tOn)^^ ^9(^(0)^(0)) (3_g4) 

and now we identify the powers of N with equation J3.58[ For any given p > 0, we must 
have: 

(m + n)! g-+^"e(nW, t^) 

= E E E ;irT^^"'' ■ ■ ■ "®<"'°'- P-''^> 

r hi • 1- • • • r • 

where in the first sum, the indices are such that 

m n 

p = 2^A:, + 2^j,-m , > 0, > (3.66) 

i=l 1=1 

and in the second sum 

r 

p = ^{2hi + h-2) , li>0, 2- 2hi - k < (3.67) 

1=1 

This equation defines m'-'^-* and t^''^ recursively in a unique way. 

Indeed, assume that we have already computed u^^\ . . . ,u 
Choose p = 2g — 1 in eq j3.65l 

r k hi ^ 



- y y ^^1^^' n^'^^) . . . w^^™) t(^-^) . . . ar'"e(w(°), tW) 

(3.68) 



m!n! 



where in the first sum we have 2g — 1 = X]I=i(2^j + h. — "2) , h > 0, 2 — 2hi — Zj < 0, and 
in the second sum we have 2g — 1 = 2 J27^i ki + 2 J27=i ji ~ which implies ki < q 
and ji < q, i.e. all the terms in the RHS are known from the recursion hypothesis. We 
have thus determined u^''\ Then, let p = 2q, we have: 
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r li hi 



- E E 



{m + n)\ 



^(fci) . . . ^{k^) tin) . . ^On) 5-+2nQ(^(0)^ ^(0)) 



m\n\ 



(3.69) 



where in the first sum we have 2g = X]I=i(2^i + h — 2), > 0, 2 — 2/;,j — Zj < 0, and 
in the second sum we have {m,n) ^ (0,1), 2g = 2^^-^/cj + 2^"^-^jj — m, which 
imphes ki < q and jj < g, i.e. all the terms in the RHS are known from the recursion 
hypothesis. We have thus determined t'^'^\ 

Therefore we have: 



It would be interesting to understand how this formula matches the tau-function 
obtained from integrability properties [2]. 

4 Holomorphic anomaly equations 

One may observe that all the terms with even powers of in formula eq J3.58l have al- 
ready appeared in another context, in topological string theory [34], and more precisely 
the so called "holomorphic anomaly equations" [8]. 

Holomorphic anomaly equations are about modular invariance versus holomorphic- 



Let us briefiy sketch the idea. String theory partition functions represent " integrals" 
over the set of all Riemann surfaces with some conformal invariant weight. In other 
words, they are integrals over moduli spaces of Riemann surfaces of given topology, and 
topological strings are integrals with a topological weight, they compute intersection 
numbers (see [36], [34] for introduction to topological strings). 

Moduli spaces can be compactified by adding their "boundaries", which corre- 
spond to degenerate Riemann surfaces (for instance when a non contractible cycle 
gets pinched). The integrals have thus boundary terms, which can be represented 
by 5-functions, and (5-functions are not holomorphic. In other words, string theory 
partition functions contain non-holomorphic terms which count degenerate Riemann 
surfaces. 

On the other hand, if one decides to integrate only on non-degenerate surfaces, one 
gets holomorphic patition functions, but not modular invariant, because the boundaries 



Z(7) = Z{t*) e NF^ + ^N^"^^u^^\nxT + ^ A^-^^t^^') 



k j 




(3.70) 



ity. 



16 



of the moduli spaces are associated to a choice of pinched cycles. Modular invariant 
means independent of a choice of cycles. 

To summarize, the holomorphic partition function is obtained after a choice of 
boundaries, i.e. a choice of a symplectic basis of non contractible cycles Ai fl Bj = 6ij, 
and cannot be modular invariant. The modular invariance is restored by adding the 
boundaries, but this breaks holomorphicity. 

There is thus a relationship between holomorphicity and modular invariance. 

Let Fg be the partition function corresponding to the moduli space of non- 
degenerate Riemann surfaces of genus g, i.e. Fg is holomorphic but not modular 
invariant (it assumes a choice of a basis of cycles Ai, Bi, i = 1, . . . ,g), and let Fg 
be the partition function including degenerate surfaces, i.e. non holomorphic but mod- 
ular invariant. The holomorphic anomaly equation discovered by 0, states that: 

dFg =\dK + |j j (4.71) 

where k is the Zamolodchikov Kahler metric symmetric matrix: 

k={t-t)-^ (4.72) 

It was found in [HI [D [I9] that: 



gEg N^-^^F, ^ ^ ^ ^ 
I k k>0,^^^^_ 

<h) pih) (r,l _ -I Nil „./ 



NT.^{'i-'ihi-l^) 



F^^>...F^^> (2/-1)!! 



(4.73) 

Remember that we use tensorial notations, and 

(2/-l)!!«:' Fi[^K . . fH^^^ (4.74) 

means in fact a sum of (21 — 1)!! terms containing all the possible pairings of 21 indices 
of the matrix k, with the 21 indices of the tensors Fj^^K 

For example to order A^~^, i.e. g = 2 we have: 

F, = F, + k(^ + ^] + 3k' (^ + 2^1^)+ 15k' (^^] (4.75) 



^2 2 J V 4! 2 3! y V2 3!3!^ 

where the last term 15 k' {Fq Y contains two topologically inequivalent types of pairings 
among the indices: 

«1,«2,«3.«4,«5,«6 L J. i 4 s c, 
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This equation can be diagrammatically represented as follows [1]: 

2 2 8 2 




(4.77) 

where each term represents a possible degeneracy of a genus 2 Riemann surface (imagine 
each link contracted to a point). The prefactor is 1/ ^Aut, i.e. the inverse of the number 
of automorphisms, for instance in the last graph we have a Z2 symmetry by exchanging 
the 2 spheres, and a symmetry from permuting the 3 endpoints of the edges, i.e. 
12 = #(Z2 X (73) automorphisms. 

Formally, eq l4.73l is very similar to eq J3.58[ with the identification: 

(2A;- 1)!! k'' 6(2^=) (4.78) 

proof: eq J4.73l is the Wick theorem expansion of the following integral fil 119]: 



dr] 



F(n)-N^7j-e*)F^-^{ri-e*)F^'iv-e*)-NH7riv-e*)K-H7j-e*) 



Z{e*) I c/r/e^'>"^'^>i-V2^^^i''(^-^*)-^'*-(''-^*)'^"'(''-^*) 

k li>0 hi>l-li/2 

dr] [t] - e*)^^' ^-NH^{n-e*)K-Hn-^') (4 



I.e. 



(4.80) 



U=0±=~-kK 



1.— 1 



where /(m, t) is nearly the same as except that we have an integral instead of a sum 
over integers: 



f{u, t) = j de e^(^-^*)" ^NH.-e*)tie-e*) ^2^nN . 
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(4.81) 



It also satisfies: 
It is clear that: 







U = 0,t=-iK-l 

which proves our claim eq l4.78[ 



dt 



9V 



{21 - 1)!! 



u=0,t= 



(4.82) 



(4.83) 



This analogy between convergent integrals obtained by summing over filling frac- 
tions, and holomorphic anomaly equations is puzzling, and it would be worth getting 
some understanding of that fact. 



5 Background independence 

So far, e* was chosen arbitrary, and eq J3.58[ eq l3.70l and the property 14.781 hold in- 
dependently of the choice of e*. Indeed 2' (7) does not depend at all on a choice of 
e*. 

If we take eq J3.58l as a definition of a string theory partition function, it seems at 
first sight that it depends on e*, but in fact it does not. Those facts are related to the 
so-called "background independence" problem in string theory [38] . 

From [6], it should be expected that if we choose e* such that the spectral curve 
has the Boutroux property: 

Boutroux property : VC , Re ^ ydx = (5.84) 

then, the formal series N^^^^Fg as well as the 6-sums in eq J3.58l and eq J3.70l should 
be convergent series, and thus we really have an asymptotic expansion instead of only 
an asymptotic series. However, this fact is not proved yet (except for the 1-matrix 
model). 

Boutroux curves in particular, are such that: 

e* = — iydxeMP , Re F^^ = Re f ydx = (5.85) 

Boutroux curves can be obtained as follows: Notice that ReF^' = —tc Imr < (the 
imaginary part of the Riemann matrix of periods is always positive, see f27\ [2S]). and 
thus — ReFo is a convex function on e* G M^, therefore it has a minimum in each cell of 
the moduli space. The minimum clearly satisfies eq J5.85[ In other words there should 
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be one Boutroux curve in each cell of the moduli space of the spectral curve. One may 
expect that each cell corresponds to a possible connectivity pattern of the generalized 
path 7. 

Notice that if the potentials are real, and the filling fraction e* is real, then Fq is 
real as well, and the Boutroux condition becomes Fq = 0. 

6 Conclusion 

In this article, we have improved the asymptotic (conjectured) formula of [9j for matrix 
integrals to all orders. We have also found an interesting connection between this 
expansion and combinatoric geometry of degenerate Riemann surfaces, through the 
holomorphic anomaly equation. 

The relationship between higher genus ^ > formal matrix integrals and nodal dis- 
crete surfaces was already known [IHIE], and here we see that there is also a relationship 
with nodal Riemann surfaces. In fact, so far all intersection numbers in Kontsevich in- 
tegral [22], or Weil-Petersson volumes [251 [26], were computed with genus zero {g = 0) 
spectral curve formal matrix models. This works tends to show that higher genus 
spectral curves have to do with nodal surfaces. This relationship needs to be further 
investigated. 
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